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Abstract. 

The monoclinic crystal structure of /3-Ga2 03 leads to significant anisotropy of the 
thermal properties. The 2a;-method is used to measure the thermal diffusivity D in 
[010] and [001] direction respectively and to determine the thermal conductivity values 
A of the [100], [010] and [001] direction from the same insulating Mg-doped 
single crystal. We detect a temperature independent anisotropy factor of both the 
thermal diffusivity and conductivity values of dl[oio]/^[ooi] = "^[oio]/^[ooi] = 1-4 ± 0.1. 
The temperature-dependence is in accord with phonon-phonon-Umklapp scattering 
processes from 300 K down to 150 K. Below 150 K point-defect-scattering lowers the 
estimated phonon-phonon-Umklapp-scattering values. 
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Semiconducting transparent oxides like /?-Ga203 provide a huge potential for future 
use in high-power electronic applications [1] and opto-electronics, like thin-hlm 
electroluminescent displays |2] or transparent held effect transistors [3]. Monoclinic jS- 
Ga203 can be grown from the melt and it is the most stable phase in a large temperature 
range from 4.2 K to 1725 K out of hve different modihcations {a - e). Recently, 
investigations of optical, electrical mn and thermal properties PE] have been carried 
out, but there are still unsolved issues. In particular, the thermal properties are yet to 
be investigated. 

Previous measurements on Gzochralski grown /d-Ga203 crystals were carried out 
for the temperature-dependent thermal conductivity in [100]-direction from room 
temperature down to 25 K, [6] and for the [010] direction from room-temperature 
up to 1200 K 1^. The room temperature thermal conductivity was found to be 
anisotropic: A[ioo] = (13 ± l)Wm“^K“^, [3] A[oio] = 21 Wm“^K“^ [7]. Additionally, 
recent investigations of the thermal conductivity of Sn-doped edge-dehned him fed 
grown (EFG) /3-Ga203 crystals were carried out for all directions from diherent samples 
using the time-domain thermorehectance method [8] with room temperature thermal 
conductivity values of A[ioo] ~ (10.9 ± 1.0)Wm“^K“^, A^oio] = (27.0 ± 2.0)Wm“^K“^ 
and A[ooi] ~ 15Wm“^K“^. A theoretical investigation from Santia et al. [9j calulated 
thermal conductivity values from hrst principles with Atheo,[ioo] = 16 Wm“^K“^, 
Atheojoio] = 22 Wm“^K“i and Atheojooi] = 21 Wm“iK“i at 300 K. 

This work provides temperature-dependent values of the thermal dihusivity in [010] 
and [001] direction as well as thermal conductivity in [100], [010] and [001] direction of 
Mg-doped insulating monoclinic /d-Ga203 bulk crystals. We implemented the 2u method 
in order to determine the thermal dihusivity along diherent crystal axes using current 
heating techniques. The complete data set was determined at the same sample. The 
method proves powerful in order to measure the thermal properties of single crystals 
along arbitrary directions. The measurements were performed in a temperature range 
from 64 K to 300 K. By applying the electrical line heater 2ci;-method we obtain 
anisotropic values from the very same crystal. Only recently, this 2ci;-method was 
presented by Ramu et al. [10] to determine the anisotrop thermal properties by using the 
transport of temperature oscillation between two parallel metal lines. Here, this method 
is hrst applied to /9-Ga203. We demonstrate it to be more powerful than the 3a;-method 
because of its accurate determination of the thermal dihusivity and conductivity for 
specihc directions. The doping with magnesium leads to electrical insulation [7|. This 
allows us to electrically separate the line heaters from the investigated crystal. From the 
temperature dependence on the thermal properties we obtain insight into the phonon 
scattering processes. 
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Figure 1. a) An arrangement of two line heater pairs on top of the Mg-doped 
/3-Ga2 03-crystal. b) The anisotropic thermal conductivity measurement setup to 
measure the temperature flow between two metal lines (Au) through the j5-Ga.20^- 
crystal with a thickness of t = 0.5 mm. The half metal line width w = 35 /rm, distance 
of the inner voltage contacts L = 1 mm and the distance between heater and sensor 
d = 240 /im are determined with optical microscopy. 

2. Materials and Method 

Our sample specific measurement setup consists of two elements: The Mg-doped 
insulating bulk / 9 -Ga 203 single crystal and two metallic heater lines, which are deposited 
on top of it. These lines can serve as heater and temperature sensors individually. 

The crystal was grown from melt by the Gzochralski method with use of an iridium 
crucible and a dynamic, self-adjusting growth atmosphere to minimize decomposition 
of Ga 203 and oxidation of the iridium crucible. The crystal is grown along the [010]- 
direction with a diameter of 22 mm or 50 mm. From the bulk crystals epi-polished 
(lOO)-oriented wafers were prepared. A detailed description of the growth process can 
be found in references mm- With an Mg-concentration of about UMg = H wt.ppm. 
the sample becomes electrically insulating. This is due to the compensation of oxygen 
vacancies (donors) with Mg [TT]. The very low concentration of Mg will only affect the 
lattice thermal properties at low temperatures. Our previous measurements [6] have 
shown, that there is no difference in the thermal conductivity of doped or undoped 
/ 3 -Ga 203 from 150 K to 300 K. 

Metal-lines were patterned by laser lithography (positive resist AZEGI 3027) and 
an Au-£lm (50 nm) was sputtered. Subsequently, a lift-off was performed with acetone 
in an ultrasonic bath. Figure a) shows two typical heater line structures. 

The metal heater lines are operated by electrical joule heating using an AG current. 
In order to determine the generated temperature oscillations AT, the temperature 
dependent resistance R = Ro + R^aAT of the metal lines is measured. Here, Rq depicts 
the resistance at the bath temperature Tq. The temperature coefficient a is determined 
by a = (1/To) ■ OR/dT. The temperature dependence of the resistance of metals is 
given by the Bloch-Griineisen-Law na 

/ T \ ^ r^o/T 5 

fl(r)-fl(4,2K)cc(^) I (1) 

As shown in figure the measured heater resistance R is in accord with equation [Tj 
From the numerical differentiation of the fitted values from equation we determine 
(y{T) in the full range of temperatures of our investigations (64 K to 300 K). 
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T[K] 


Figure 2. The temperature dependent resistance of the four heater lines approximated 
with equation]^ the so called Bloch-Griineisen-law. The temperature coefficient a{T) 
for both heater lines is calculated from the numerical derivation of the fit curve. 


Joule heating is achieved with an alternating power P = U-I with the oscillating voltage 
U = Uq ■ cos{ut) and an oscillating current I = Iq ■ cos{ut). This leads to an oscillating 
temperature change AT oc P = {UqIq/2) ■ (1 + cos{2ut)) (see ref[T3]) with the angular 
frequency uj = 2nf and / denotes the measurement frequency. These temperature 
oscillations are transported within the crystal according to its thermal properties. To 
detect temperature oscillations in a metallic sensor line in a distance d parallel to 
the heater with a direct current /q.dc, a resistance Ps and a temperature coefficient 
as, one needs to detect the correlated resistance oscillations R^{2uj) oc aATs(2a;) 
and therefore the measureable voltage oscillations of two times the heater frequency 
U2to{2uj) = Is,DC • Ps(2a;). Here, voltage oscillations are detected as a function of the 
transported temperature oscillations ATs by 


U2.= 


as-fs,DcPsV^ 


ATs 


( 2 ) 


A Lock-In ampliher (SR830) is used to detect these higher harmonic voltages. A 
schematic and description is given in hgure 
From the frequency dependence of the voltage signal P 21 V the thermal properties can 
be determined as described below. 


3. Results 

The determination of the thermal diffusivity and conductivity by the 2a;-method is 
depicted as follows. The thermal conductivity A, thermal diffusivity D and specihc heat 
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Figure 3. The 2a;-method is based on two seperate electrical circuits, one heaterline 
with dominating alternating heating current and the sensorline with direct current and 
heat induced second harmonic voltage U 2 u:- 


capacity Cy are related to each other by the equation 
X = D-Cv . 


(3) 


The thermal conduction differential equation 
d^AT{r,t) 1 dAT{r,t) _ 
dr'^ ^ D dt 

can be solved [13] for an oscillating temperature AT to 


(4) 


AT(r) = —-- ■ K,{qr) (5) 

^ Theater^ 

with P/Theater the power per unit length of the heater, Kq the zero-order Bessel function 
of the second kind and q represents the frequency dependent inverse thermal penetration 
depth q{f) = ^jAirf/D. In order to determine the thermal diffusivity D along the crystal 
axes [010] and [001], equation has to be utilized to determine the heat flow between 
the heater and the sensor lines by averaging their widths |T5| 

P 1 r'^heater ^ P'^lsensor 

^ - "" i^o((g[x]-(d + o-p))dodp 

z,u/heater d-tuheater ^ “^sensor '^^sensor 


Here, we denote the inverse thermal penetration depth as q[x\ = the half 

heater line width as tCheater, the half sensor line width as Wsensor, the heating power as P, 
the length of the heater line as Theater and the average distance between the heater and 
the sensor line as d. The integration parameters o and p depict the distance along the 
line heater widths. The effective thermal conductivity is A = y/Ap] • A[y] (see ref. [T0]h 
An exemplary measurement of AT(/) oc U 2 ^{f) is shown in figuretogether with a fit 
using equation by varying the thermal diffusivity D[x] and the thermal conductivity 
A. The integrals were solved numerically. The axis [x] denotes the direction from the 
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heater to the sensor line and [y] the direction perpendicular to the surface plane, here 

[ 100 ], 

Here, two heater and sensor lines are placed along the [001] and [010] direction. 
Therefore, the thermal diffusivity values /l[oio] and Il[ooi] are directly obtained (hgure 
5), as well as the effective thermal conductivity values -^Apoo] • A[oio] and -^Aiioo] ■ A[ooi] 
shown in hgure 

With the knowledge of the isotropic specihc heat capacity Cy one can determine 
the thermal conductivity values for all three primary axes as shown in table 



f [Hz] 


Figure 4. The 2a;-voltage oscillations were measured with displaced heater and sensor 
line perpendicular to the [001]-direction. The approximation was performed using 
equation!^ The plot was by way of example taken at 250 K. 


4. Discussion 

From the 2a;-measurements we obtained directly the temperature dependent thermal 
diffusivities in [010]- and [001]-direction as shown in hgure We observe higher thermal 
dihusivity values in [010]-direction then in [010]-direction. For high temperatures from 
150 K to 300 K we observe a linear dependence in logarithmic scale for both directions 
with D{T) = a ■ T™ (solid lines). Both linear approximations have the same slope 
of m = 1.9 ± 0.1 which means the anisotropy factor is D[oio]/D[ooi] = 1-4 ± 0.1 and 
temperature independent. The temperature dependence of the dihusivity is determined 
by the temperature dependence of the phonon mean free path Aph and the temperature 
independent velocity of sound Vg as 

D(T) = ^Aph(T)-n, 


(7) 
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Figure 5. Temperature-dependent thermal diffusivity of insulating / 3 -Ga 203 on a 
logarithmic plot vs. the inverse temperature. The approximaton D(T) = a ■ T™ (solid 
line) reveals the same slope m for both directions. The measurements agree with 
phonon-phonon-Umklapp scattering (equation]^ dotted line) in the high temperature 
regime from 150 K to 300 K. Below 150 K, point defect scattering presumably lowers 
the diffusivity values. 


The temperature dependence of the phonon mean free path on behalf of phonon-phonon 
interaction with Umklapp-scattering would be [16] 

A oc - 1) (8) 

with as Debye-temperature. The approximation of the diffusivity values with 
equation (dotted lines) fits well in the temperature range from 150 K to 300 K. Below 
150 K point-defect-scattering lowers the estimated thermal diffusivity presumably by 
inhibiting the phonon mean free path with an additional scattering effect. The Mg- 
concentration dominates all other defects and is most likely the cause of point-defect- 
scattering. 

From the experiment we determine the temperature dependent thermal 
conductivity effective values Aac = ■\/A[ioo] • A[ooi] and Aab = \/A[ioo] • -^[oio] as shown in 
hgure The anisotropy factor of the two effective thermal conductivity values squared 
should be the same as the anisotopy factor of the thermal diffusivities because 

\/A[ioo] ■ A [ 010 ] 

\/-^[ioo] ■ -^[ooi] 

The temperature depedence of Aab and Aac can be seen in hgure approximated by 
the equation A(T) = b-T'^ with the slope n = 1.3±0.1. For high temperatures T > 200 K 
this slope value n is found both experimentally by Guo et al. [8] and theoretically by 



= lm = i,4±0.1 (9) 

A [001] 
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Figure 6. Temperature dependent thermal conductivity of insulating /?-Ga 203 on a 
logarithmic plot with a A(T) = b ■ T” approximation. 


axis 

L.p. dn 

D 

A 

Aex,ref [H] 

Atheo,ref [9] 


A 

mm^s“^ 

Wm-^K-i 

Wm-^K-^ 

Wm-^K-^ 

a [100] 

12.2 

3.7 ±0.4 

11 ± 1 

10.9 ± 1.0 

16 

b [010] 

3.0 

9.6 ±0.5 

29 ±2 

27.0 ±2.0 

22 

c [001] 

5.8 

7.1 ±0.4 

21 ±2 

15 

21 


Table 1. Lattice parameter (L.P.), thermal diffusivity D and thermal conductivity A 
values for the different crystal axis of our Czochralski grown j3-G&2^3 and literature 
values Aref at room temperature. 


Sanita et al. [9]. However, for low temperatures the dominance of phonon mean free 
path limiting scattering processes is dependent on the crystal quality and differs. For 
room temperature using the specific heat Cy = 0.5 kJkg“^K“^ = 3 kgcm“^K“^ one 
can calculate all main axis thermal conductivity and diffusivity values as listed in table 
The crystallographic order of the thermal conductivity values is in agreement with 
experimental work by Guo et al. as well as the theoretical result of Sanita et al. [9]. The 
lowest thermal conductivity value is along the [100] axis, the highest for the [010] axis. 
The values of A[ioo] and A[oio] are in accordance with Guo et al. [8] the value of A[ooi] 
is in our case about 30% higher than in that of Guo et a/. [8]. The anisotropy factors 
including the [100]-direction are A[oio]/A[ioo] = 2.5 ± 0.3 and A[ooi]/A[ioo] = 1-9 ± 0.2. 
Additionally performed 3a; measurements conhrm the presented values with a lower ac¬ 
curacy and are shown in the Appendix. 
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The temperature-dependent values of the thermal diffusivity in [010] and [001] direction 
as well as thermal conductivity in [100], [010] and [001] direction of Mg-doped insulating 
monoclinic / 3 -Ga 203 bulk crystals were measured by the 2a;-method. We hnd an 
anisotropic behavior for the thermal diffusivity and conductivity along [100], [010] and 
[001]. The anisotropy factor T>[oio]/^[ooi] = A[oio]/A[ooi] = 1.4 ± 0.1, A[oio]/A[ioo] = 
2.5 ± 0.3 and A[ooi]/A[ioo] = 1.9 ± 0.2 are temperature-independent. The room 
temperature values for the thermal conductivity in the main crystal axes are determined 
to A[ioo] = (11±1) Wm“^K“^, A[oio] = (29±2) Wm“^K“^ and A[ooi] = (21±2) Wm“^K“^. 
The temperature dependence of the thermal conductivity and diffusivity conhrms the 
model of phonon-phonon Umklapp scattering for T > 150 K. At lower temperatures 
the estimated thermal diffusivity and conductivity values are reduced by point-defect 
scattering. 
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Appendix 


Usually 3a; measurements are used to determine the thermal conductivity of materials. 
To conhrm the presented values for the thermal conductivity one can compare the 2u- 
method used in this article with the 3c;-method|13j measurement results. Here, we 
present measurements on a similar Mg-doped insulating / 3 -Ga 203 single crystal with 
two heater lines to determine the effective thermal conductivity y^A[ioo] x A[ooi] and 
compare it with the values of the 2a;-method shown in hgure 

To determine the effective thermal conductivity with a single heater line one 
has to detect the temperature oscillations AT{2u) induced by the heater lines AG 
power P{2u). These temperature oscillations lead to resistance oscillations in the 
heater line R = Rq + RoaAT{2(jj). The measured voltage U contains two parts 
U = R{2oj) ■ /(la;) = U{loj) + U{3uj) with 


U{3u:) = 


( 10 ) 


To detect the rather small voltage U{3u) within a huge voltage U{lu) a Lock-In ampliher 
and a hlter were used. The frequency dependence of the temperature oscillations is 

sin(/tn) 


p r 


kw\Jk? -|- q- 


-Ak 


( 11 ) 
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Figure 7 . Exemplary plot for the frequency-dependent voltage fitted with the 

approximation of equation 11 solid lines). 


Here, we denote the inverse thermal penetration depth q = the half heater 

line width ta, the heating power P, the length of the heater line Lheater- Here, the 
effective thermal condnctivity[T^ is A = y/Apoo] ■ A[ooi]. 

The agreement of the measnrement with this solntion is shown in hgure The 
direct comparison of the obtained values of the 2aj and the So; measurements are shown 
in hgure 

There is a good agreement between these two measurements with a better accuracy 
of the 2u: measurement. The slopes n of the approximation A = h-T'^ are n 2 u) = 1.3±0.1, 
= 1.4 ±0.1 and = 1.3 ±0.1 are the same. The room temperature (T = 297K) 
values are A 2 a; = 15 ± 1, As^q = 14 ± 2 and = 13 ± 2. There are three reasons why 
the 2u: method is prefered to the 3n; method. First, we get well-dehned thermal diffu- 
sivity values for specihc directions. Second, there are no parasitic 2u: signals within the 
current, but every small heat source produces a den signal. Third, the Lock-In ampliher 
does not have to take care of an oscillating Ui^i signal. 
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